Landau model for the elastic properties of the ferroelastic crystal Rb 4 LiH 3 ( 804)4. 
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We report a detailed investigation of the elastic properties of Rb4LiH3(SC>4)4 measured as a 
function of temperature and pressure. Rb4LiH3(S04)4 is known to show a 4 — > 2 ferroelastic phase 
transition at T c = 134 K. In order to clarify the nature of the order parameter associated with 
that structural transition, we compare our finding to two distinct phenomelogical Landau models. 
The coupling parameters of both models are all determined using the temperature dependence of 
the strain tensor, as well as the pressure dependence dT c /dP = 19.1 ± 0.2 K/kbar, prior to the 
calculation of the elastic constants. Our comparison indicates that the ferroelastic transition in 
Rb4LiHs(S04)4 is fully consistent (within a few percent) with the predictions of a pseudo-proper 
model, showing at the same time that the primary mechanism leading to the phase transition is 
not driven by strains. Our analysis also refutes the idea that Rb4LiR"3(S04)4 is the first ferroelastic 
compound showing incomplete softening, of the effective soft acoustic mode C' so f t , at T c . 

PACS numbers: 



I. INTRODUCTION 



Over the years, there has been a series of con- 
troversies regarding the properties of the ferroelas- 
tic compound Rb 4 LiH3 (804)4. Initially designated as 
LiRb 5 (S04)2l.5H 2 S0 4 1 i 2 i 3 i 4 , subsequent chemical analy- 
sis revealed that Rb4LiH 3 (804)4 is the proper descrip- 
tion. It was also suggested that Rb4LiH3 (804)4 belongs 
to the 4mm point group in the paraelastic phase at ambi- 
ent temperature. Consequently, a 4mm — > 2mm symme- 
try change, associated with a ferroelastic phase transition 
at 134 K, was assumed^. Additional X-raj^ and neu- 
tron diffraction 8 measurements finally established that 
the proper group-subgroup symmetry change is actually 
4^2. Based on this new phase sequence, the elas- 
tic properties of Rb4LiH 3 (S04)4 have been revisited us- 
ing Brillouin scattering^ and ultrasonic^ measurements. 



Both experimental approaches now confirmed that the ef- 
fective elastic constant (Cn — Ci2)/2 softens as the tem- 
perature is reduced to T c , which is perfectly compatible 
with a 4 — > 2 phase transition. Nonetheless, in order 
to explain their respective findings, each group proposed 
models which are fundamentally different. On the one 
hand, the ultrasonic measurements^, which show a sud- 
den drop at T c in the velocity of longitudinal modes prop- 
agating along [100] and [001], are interpreted within the 
framework of a pseudo-proper ferroelastic model. Con- 
sequently, the authors assume that the order parame- 
ter is not one of the strains or a combination of strain 
components. However, in that case the order parameter 
must at least exhibit the same symmetry as the strain 
(ei — €2) or egii. On the other hand, Mroz et ali^ claim 
that Rb4LiH3 (804)4 shows incomplete softening at T c 
and proposed a model which they claim is consistent 



with that observation. In their paper, Rb4LiH3 (804)4 
is described as a proper ferroelastic compound where the 
strain combination e so ft = ot\ (ei — 62) + ae ee acts 
as the order parameter. Incomplete softening in ferroe- 
lastic crystals is not a common phenomena, but it has 
previously been observed in BiVO^pi^. In that case, the 
softening is mediated by the coupling between optical 
and acoustical modes. Considering that the model pro- 
posed by Mroz et al£ includes no such coupling term, 
we believe that their conclusion regarding the soft mode 
Rb4LiH3 (804)4 requires further investigation. Thus, in 
order to clarify the nature of the ferroelastic transition, as 
well as, to validate or refute the possibility of incomplete 
softening in RbiLiH3(S04)4, we present a detailed inves- 
tigation of the elastic properties of Rb4LiH3(S04)4 as a 
function of temperature and pressure. We also compare 
the predictions of a proper and a pseudo-proper ferroe- 
lastic models in order to ascertain the character of the 
ferroelastic transition in Rb^LiHs (804)4. 



II. EXPERIMENT 

The Rb4LiH3(S04)4 crystal was grown by the Crystal 
Physics Laboratory of the Institute of Physics at Mick- 
iewicz University, Poland^. For the ultrasonic investi- 
gation, several samples in the form of cubes of about 
3x3x3 mm 3 were used to measure the sound velocity 
along different crystallographic directions. Longitudinal 
and transversal waves were generated and detected using 
30 MHz lithium niobate transducers mounted on one face 
of the crystal. The relative variations in the sound veloc- 
ity (AV/V) were obtained using a high-resolution pulsed 
acoustic interferometer as a function of temperature and 
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FIG. 1: Temperature dependence of the sound velocity of 
longitudinal (L) and transversal (T) modes measured along 
the principal crystallographic directions of Rb4LiH3(S04)4. 
The experimental results are represented by open symbols 
while the continuous lines correspond to predictions derived 
from the pseudo-proper ferroelastic model presented in the 
Appendix. 



Fig. n and Fig. [21 our high resolution sound velocity re- 
sults obtained for RlxiLiHa (804)4 as a function of tem- 
perature and pressure using modes propagating along 
the principal axes. In these figures, the experimental 
results are represented by open symbols while the con- 
tinuous lines correspond to predictions derived from the 
pseudo-proper ferroelastic model presented in the Ap- 
pendix IVII Bl Let first point out that our results are 
very consistent with previous studie s 9 ' 10 realized as a 
function of temperature. As for the Brillouin scatter- 
ing measurements^, the results presented in Fig. ^ in- 
dicate that the largest softening is observed for trans- 
verse modes propagating along [110] with a [TlO] polar- 
ization (V r p- 10 ] [110]). In addition, our high resolution 
measurements show how the other modes change in the 
vicinity of the ferroelastic transition (T c = 132.8 K). 
In particularly, we note that the velocity of longitudi- 
nal modes presented in Fig. ^ drop by a few percent at 
T c . This observation is also consistent with previous ul- 
trasonic measurements^. Our investigation is comple- 
mented with the first measurements realized as a func- 
tion of pressure which we presented in Fig. |2J At ambi- 
ent temperature, the data indicate a phase transition at 
a pressure of P c = 8.6 ± 0.2 kbar. Considering that the 
observed anomalies on the elastic modes at T c and P c are 
very similar, we can assume that the transition observed 
at P c = 8.6 ±0.2 kbar corresponds to a 4 — > 2 ferroelastic 
structural transformation. To determine the pressure de- 
pendence of the critical temperature, we have carried out 
a series of sound velocity measurements as a function of 
temperature at different pressures. The results presented 
in Fig. |31 have been obtained using longitudinal waves 
propagating along the z-axis. From these measurements, 
we find that the ferroelastic transition temperature in- 
creases at a rate of dT c /dP = 19.1 ± 0.2 K/kbar with 
pressure (see inset of Fig. [2Jl - 



pressure. For measurements realized at high pressure, the 
transducer-sample assemblage was inserted in a Cu-Be 
pressure cell filled with a 3-Methyl-l-butanol fluid acting 
as the pressure-transmitting medium. A small wire of 
Lead mounted close to the sample was used during the 
experiments to determine the actual pressure at different 
temperatures. 



III. EXPERIMENTAL RESULTS 

In this paper, the labels L and T are used to identify 
longitudinal and transverse modes, respectively. In addi- 
tion, the Miller index adjacent to the label T gives the po- 
larization of the mode while the other index corresponds 
to the direction of propagation. For example, Vr[ooi] [100] 
stands for the velocity of transverse wave propagating 
along the crystallographic direction [100] with its polar- 
ization along [001]. Using that notation, we present in 



IV. DATA ANALYSIS 

The results presented in Fig. |3] show that the tem- 
perature dependence just below the critical temperature 
changes significantly between ambient pressure and a 
moderate pressure of P = 1.5 kbar. We believe that 
this might be an indication that the observed temper- 
ature dependence in the ferroelastic state is not purely 
intrinsic. One of the characteristic of ferroelastic com- 
pounds is that structural domains appear in the ferroe- 
lastic phase. Moreover for soft ferroelastic materials, it is 
possible to switch the orientation of these domains by ap- 
plying an uniaxial stress. In the case of Rb4LiH3 (804)4, 
the domains have been observed^ and consist of two mu- 
tually perpendicular walls. Considering that the pressure 
is not perfectly isotropic in pressure cell, the unusual be- 
havior of the velocity close to T c might be associated 
with a modification of the domain structure with increas- 
ing pressure. Consequently, before attempting to analyze 
any results in the ferroelastic phase, it is crucial to de- 
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FIG. 2: Pressure dependence of the sound velocity of lon- 
gitudinal (L) and transversal (T) modes measured at room 
temperature along the principal crystallographic directions of 
Rb4LiH3(SC>4)4. The experimental results are represented by 
open symbols while the continuous lines correspond to pre- 
dictions derived from the pseudo-proper ferroelastic model 
presented in the Appendix. 



tcrmine what might be the influences of these domains 
on the measured quantities. Thus, in order to ascertain 
if our measurements represent the intrinsic elastic prop- 
erties of Rb 4 LiH 3 (S04)4, we compare in Fig. ^the tem- 
perature dependence of the velocity of longitudinal waves 
propagating along two orthogonal directions in the ab 
plane ([100] and [010] directions). As the x and y direc- 
tions are no longer equivalent in a monoclinic phase, one 
would expect a significant difference between these two 
modes. The fact that both sets of data are practically 
identical below T c indicates that the size of the domains 
is smaller then our acoustic wavelength (A « 100 \irn) at 
30 MHz. Thus, based on this comparison, it is clear that 
our results obtained in the ferroelastic phase, as well as 
those obtained by Breczewski et aliiS, do not reflect the 
intrinsic elastic properties of Rb 4 LiH3(S0 4 ) 4 . For that 
reason, the Landau models presented in the Appendix are 
tested using results obtained principally in the paraelas- 
tic phase (tetragonal phase). 
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FIG. 3: Temperature dependence of the sound velocity of 
longitudinal waves propagating along the z-axis (Vl[001]) 
measured at different pressures. The inset shows the pres- 
sure dependence of the ferroelastic critical temperature T c for 
Rb 4 LiH 3 (S0 4 )4. 

3500 



CO 

E 



3450 
3400 
> 3350 
3300 
3250 



oV L [100] 




□ V L [010] 


a 




o 
a 

o 

□ 


□ 


1 "V o 

n O 

n O 
n O 

: ° °n o 


o 

& 

o 




§ 


%; 





80 



100 



120 
T(K) 



140 



160 



FIG. 4: Comparison of the temperature dependence of the ve- 
locity of longitudinal waves propagating along two orthogonal 
directions in the ab plane ([100] and [010] directions). 



So far, the elastic properties of Rb 4 LiH 3 (S0 4 )4 have 
been qualitatively analyzed by different groups&iS using 
two different Landau-type models. In one case, they con- 
sider that the order parameter corresponds to the strains 
ei — e2 and (proper ferroelastic) while in the other sce- 
nario the order parameter is unknown but has the same 
symmetry as the strains mentioned previously (pseudo- 
ferroelastic) . In this paper, we present (see Appendix) a 
comprehensive derivation of a proper and pseudo-proper 
models for a 4 — > 2 ferroelastic phase transition. Both 
of these models are derived using a minimum number 
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of coupling factors. Our principal goal is to see if one 
could determine the true nature of the order parameter 
based on a detailed analysis of the elastic properties of 
Rb 4 LiH 3 (S0 4 )4- 

According to both models derived in the Appendix (see 
Eg. 1251 and Eq. |22Jl, the elastic constants Cn, C\e, and 
Cqq are expected to soften while C12 becomes stiffer as T c 
is approached from above (or approach P c from below). 
The main qualitative difference between these models is 
that the softening and stiffening are expected to be linear 
for a proper ferroelastic transition while it could be non- 
linear in the other case. Using the relations given in refer- 
ence [10] and our high resolution velocity measurements, 
it is easy to determine the temperature and pressure de- 
pendence of the principal elastic constants. Considering 
that the value of C\q is an order of magnitude smaller 
than Cn^iS, we can consider, as a reasonable approxi- 
mation, that the temperature and pressure dependence 
of C\\ is captured by that of Vl[100]. Thus, our high 
resolution measurements indicate that C\\ indeed soft- 
ens, however, its temperature and pressure dependence 
are significantly non-linear. This first constatation is re- 
inforced by a comparison of the velocity measurements 
performed on Rb^LiH^SO^)^ and the isomorphic com- 
pound K^LiH^SOi)^ As KiLiH^SOi)^ shows no fer- 
roelastic transition^, it can be used to accurately deter- 
mined the anharmonic temperature dependence of the 
elastic constants in the absence of softening due to the 
transition. As an example, we compare in Fig. [S] the re- 
sults obtained for longitudinal modes propagating along 
the z and x directions. Along the z-direction, we note 
that both compounds have the same absolute velocity at 
room temperature and exhibit the same temperature de- 
pendence. This is consistent with both models as no 
influence of the soft mode is expected on C33 in the 
para-elastic phase. Along the x-direction, the compar- 
ison of Vt[100] for these two compounds indicates that 
the contribution associated with the soft mode is indeed 
non-linear and thus more compatible with the pseudo- 
ferroelastic description. 

So far, our qualitative observations seem to indi- 
cate that the elastic properties of Rb 4 LiH 3 (804)4 are 
more consistent with the predictions based on pseudo- 
ferroelastic model. A more rigorous test can be per- 
formed by carrying on a detailed numerical compari- 
son using the prediction giving in the Appendix. While 
the elastic constants are determined by sound velocity 
measurements realized at room temperature, the magni- 
tude of the coupling terms, as well as a and A4, need 
to be determined independently. Moreover, it is also 
important to determine the anharmonic contributions 
which arc independent of the transition. In the case 
of Rb4LiH 3 (S04)4, these contributions have been deter- 
mined from velocity measurements performed on the iso- 
morphic compound K^LiH^SOij^. Finally, let men- 
tion that both models include only eight adjustable pa- 
rameters and that six of them are determined indepen- 
dently using the temperature dependence of the strains, 



the pressure dependence of T c (dT c /dP = 19.1 K/kbar), 
and the normalization of the order parameter. The last 
two parameters £ and 77, which only influence C44 and 
C45 , can be set manually. To illustrate how stringent this 
process is, we present in Fig. the predictions based on 
Eas.l9"l -inland the temperature dependence of the strains 
obtained from high-resolution neutron diffraction results 
published by Mroz et al£ As ei — e-i and show a well de- 
fined mean field temperature dependence, we can safely 
assume that the neutron data reflet the microscopic prop- 
erties within the macroscopic domains in the ferroelastic 
phase. When all adjustable parameters are determined, 
we then calculate the velocities along different directions 
by solving the Christoffel's equation using the elastic ten- 
sors Eg. 1231 and Eq. [221 for the pseudo-proper and proper 
ferroelastic transition, respectively. The results of the 
calculations (continuous line) based on the pseudo-proper 
model are presented in Fig. lll2l along with the experimen- 
tal data points. In the para-elastic phase, the agreement 
between the prediction of the pseudo-proper model and 
experimental results, as a function of temperature and 
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FIG. 5: Comparison of the elastic properties of 
RbiLiHsiSOiji and the isomorphic compound 
KiLiHz^SOi)^ using the temperature dependence of 
Vi[100] and Vl[001]. 
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FIG. 6: Temperature dependence of the spontaneous strains. 
The points correspond to results obtained from Mroz et ali- 
while the continuous line are calculated using Eg. IS1- 1TH 

pressure, is quite remarkable (within of percent). Nat- 
urally, the agreement in the ferroelastic phase is not as 
good considering that the elastic properties are signif- 
icantly affected by the existence of structural domains 
below T c . As an ultimate test, we compare in Fig. [7| 
the predictions of the proper model (dash line) and the 
pseudo-proper model (continuous line) relative to exper- 
imental data showing the largest variation. These varia- 
tions are obtained using transverse modes propagating in 
the xy plane at an angle <fr relative to the [110] direction 
(with <j) — +5°, 0°, — 10°) with its polarization in the 
xy plane. From that comparison, it is clear that we sys- 
tematically obtain a better agreement with the pseudo- 
proper model rather than the proper ferroelastic model. 
In Fig. 0we have limited our comparison to V^r^o] [110], 
however, we reach the same conclusion using results ob- 
tained for Vl[1Q0]. 

V. SOFT MODE 

The data presented in Fig.[7|also confirms a larger soft- 
ening as the direction of propagation, for transverse mode 
polarized in the xy plane, is a few degrees away from the 
[110] direction (approximately —10° relative to the [110] 



direction). This observation complies with previous Bril- 
louin scattering measurements^ which show less softening 
as the direction of propagation is changed toward [010] 
instead of [100]. However, in their data analysis, they 
also claim that Rb4LiH3(S04)4 shows incomplete soften- 
ing at T c , a claim that we revisit in this section. The 
effective modulus for the soft acoustic mode of a 4 — > 2 
ferroelastic phase transition can be obtained by finding 
the eigenvalues of the submatrix corresponding to the 
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FIG. 7: Temperature dependence of the velocity of transverse 
waves polarized in the xy plane and propagating along differ- 
ent directions in the xy plane (0 = +5°, 0°, and —10° relative 
to the [110] direction). Here, the continuous and dash lines 
correspond to predictions based on a pseudo and proper fer- 
roelastic models presented in the Appendix. 

strain components e\ — ei and e§, 

/ C11-C12 Q \ 

\ Ci6 Cqq J 

We would like to point out that the submatrix used here 
differs from the one published by Boccarai 3 -. Consider- 
ing that elastic tensors must be symmetric, it is obvious 
that some of the results found in Boccarai 3 - are inexact 
and that they correspond to simple topographic errors. 
Thus, contrary to the expression used by Mroz et aL 9 , 
which was derived from expressions found in Boccara 13 , 
the appropriate relation for the soft modulus is 



1 I Cn - C12 ., 
t sof t - 2 2 + 1 



'C n -C 



12 



66 



- Cfi 



(2) 



while the associated direction of propagation^ is given where the angle <f> is defined relative the [110] direction, 
by Using the measured elastic constants of Rb4LiH 3 (804)4 

4Cir 
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at room temperature, the expected direction for the soft 
mode should correspond to <f> = —13°. This prediction 
(45° + 4> — 32° relative to the x-axis) agrees well with the 
orientations of the domain walls observed in the ferroe- 
lastic phased. As that direction does not coincide with 
any of the principal crystallographic directions, this has 
contributed to delay the experimental determination of 
the actual soft mode in Rb4Lifi3 (804)4. According to 
our recent investigation, the data obtained at <\> = —10° 
indicate a larger variation relative to the results obtained 
at <fi — 0°, which is very consistent with the numerical 
prediction. Unfortunately, due to a rapid increases in 
the acoustic attenuation close to T c , we were not able 
to perform the measurement down to T c . Nevertheless, 
considering that the theoretical prediction based on the 
pseudo ferroelastic model agree well with our experimen- 
tal data, we believe that within the experimental uncer- 
tainties there is no strong evidence of incomplete soft- 
ening in Rb4LiH3(S04)4. Thus, it is clear, based on 
the results presented in this section, that the conclusion 
reached by Mroz et al£ is inexact. 



in the presence of linear coupling between acoustical and 
optical modest, while the actual transition takes place 
at T c , the extrapolated temperature dependence of the 
soft optical energy in the paraelastic phase is expected to 
vanish at a temperature T D . Based on our numerical pre- 
dictions, T D is approximately equal to 20 K in the case of 
Rb4LiH 3 (804)4. We believe that further investigations 
are certainly necessary in order to clarify the nature of 
the actual mode responsible for the ferroelastic transition 
in Rb 4 LiH 3 (S0 4 )4. 

Finally, as we have clearly demonstrated in this pa- 
per, there is no strong experimental evidence that 
Rb4LiH3 (804)4 show incomplete softening of the soft 
elastic constant at T c . Consequently, the ferroelastic 
character of Rb4Lifi3(S04)4 is not necessary unique, for 
example, its properties seem very similar to those of 
(NH4)4LiH3(S0 4 )4 which shows a 4 — > 2 ferroelastic 
phase transition at 236 K 16 . It would be interesting to 
use a similar numerical analysis in order to determine the 
nature of that ferroelastic transition. 



VI. CONCLUSIONS 

In this investigation, we have presented a detailed anal- 
ysis of the elastic properties of Rb 4 LiH 3 (804)4. A sub- 
stantial set of experimental data, obtained as a function 
of temperature and pressure, is compared to two Lan- 
dau models associated with a 4 — > 2 ferroelastic phase 
transition. The possible scenarios considered here cor- 
respond to the proper ferroelastic case where the in- 
stability is driven by the strains while in the pseudo- 
proper situation the order parameter is unknown. In 
the later case, even the order parameter is unknown, its 
symmetry must be identical to the order parameter de- 
fined for the proper ferroelastic case. Both models are 
derived using a minimum number of coupling parame- 
ters which are adjusted from thermal expansion data2> 
and the pressure dependence of the critical temperature 
dT c /dP = 19.1 ± 0.2 K/kbar obtained in this work. Our 
analysis clearly shows that within the frame work of the 
pseudo-proper ferroelastic model, it is possible to obtain 
predictions that are simultaneously consistent with the 
temperature and pressure dependence of the elastic con- 
stants, the pressure dependence of T c and thermal ex- 
pansion. No such quantitative agreement is obtained us- 
ing the proper ferroelastic model. Thus, as suggested 
by our investigation, Rb 4 LiH3(S0 4 )4 should be consider 
as a pseudo-proper ferroelastic compound. If our con- 
clusion is accurate, this implies that the nature of the 
order parameter, driving the ferroelastic transition in 
Rb4LiH3(S04)4, is still unknown. So far, there has been 
a limited number of Raman scattering measurements on 
Rb4LiH3(S04)4 15 . However, even if they observed a soft 
optical B mode, the authors of Ref. 15 claim that the 
variation of that particular mode is insufficient to ac- 
count for the transition. It is worthwhile to note that, 



VII. APPENDIX 

A. Landau Models 

The direct group-subgroup relationship for the 4^2 
ferroelastic phase transition imposes that the sponta- 
neous strainii associated with this transition must cor- 
respond to the combination ot\ (e\ — e2) + ae eg. This 
effective strain can then act as the principal order pa- 
rameter (proper ferroelastic transition) or couple bilin- 
early to the order parameter Q. In the later case, the 
order parameter Q could correspond to the softening of 
an optical phonon and the phase transformation is then 
referred to as a pseudoproper ferroelastic transition 11 . 
Even if there is no direct evidence of a soft optical mode 
in Rb4LiH 3 (S04)4^, this later scenario should not be 
dismissed. As we show in this paper, these two types of 
ferroelastic transition should in principle produce subtle 
differences in the elastic properties of Rb4LiH3(S04)4. 
The predictions, derived from the models presented in 
the next section, are then compared to our experimen- 
tal data in order to determine whether Rb4LiH3 (804)4 
should be regarded as a proper or pseudo-proper ferroe- 
lastic compound. 



B. Pseudo-proper Ferroelastic Model 

Considering that we are particularly interested in the 
pressure and temperature dependence of the elastic prop- 
erties, we express the Gibbs free energy into four distinct 
contributions such as 

G(Q,e i )=F L (Q)+F el (e i )+F c (Q,e i )+F P (P,e i ) . (4) 
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Here, Fl(Q) is the usual second-order Landau- type free 
energy in terms of an order parameter Q, 



Fl(Q) 



1 



A 2 Q 



1 



A 4 Q 4 



(5) 



where A 2 — a (T—T ) depends explicitly on temperature 
while the pressure contribution is introduced in Fp(P, ei) 
of Eq. In the presence of external stresses Cj, this 
contribution is normally written as — a% ej. In the 
case of an hydrostatic pressure P the stress associated 
with longitudinal strains (i = 1, 2, 3) corresponds to 
cr, = —P while for shear strains <jj = (j — 4, 5, 6), 
hence 



F P (P, ei ) =P (ei+e 2 +e 3 ) 



(6) 



This contribution is crucial in order to calculate the pres- 
sure dependence of the elastic constants as well as to 
determine the pressure dependence of the critical tem- 
perature. Furthermore, this is achieved without adding 
any additional adjustable parameters. The elastic energy 



F e i(ei) is also straightforward to derive as it is imposed 
by the symmetry of the high temperature phase. Thus, 
for a point group 4 symmetry^ 



Fel(ei) = \ C n (e? + e%) + ~ C 33 e\ + \ C 44 {e\ + e§) 
1 



+ | C 66 eg + C12 ei e 2 + C 13 (e x + e 2 ) e 3 
+ C' 16 (ei - e 2 ) e 6 



(7) 



where Cjfe represent the bare elastic constants at high 
temperatures. Finally, the most critical term corresponds 
to the coupling energy F C (Q, ei) which takes into account 
the coupling between the strain components ei and the 
order parameter Q. For a 4 — > 2 pseudo-proper ferroelas- 
tic transition, we know that the order parameter Q must 
transform as the spontaneous strain ai(e± — e%) + ct§e§ 
under the symmetry operations of the high temperature 
phase. Thus, using fundamental symmetry arguments, 
the coupling terms considered in this paper are 



J 



F C (Q, = 0Q (ex - e 2 ) + 7 Q e 6 + 6 Q 2 e 3 + A Q 2 (ei + e 2 ) + C e 4 e 5 Q + 77 (e\ + el) Q 2 



r 



(8) 



Here, we only consider terms which are fully compati- 
ble with a 4 — > 2 structural transformation. The Gibbs 
energy, as defined in this paper (Eq. I4l8p . gives a phe- 
nomenological framework which can now be used in order 
to derive the elastic properties of a 4 — > 2 pseudo-proper 
ferroelastic compound. From the minimization of the 
Gibbs energy with respect to e,, we obtain a series of ex- 
pressions for the spontaneous strains ei(Q). At ambient 
pressure these relations correspond to 



ei - e 2 



ei + e 2 



2 ( 7 Ci 6 - g C 66 ) 

2 C 16 



(Cn — C12) 6*66 

2 (<5 C13 ~ A C 33 ) 

2 C 13 



63 = 

e 4 = 

es = 

ee = 



(Cn + C 33 
(5 (C n + Cia) 



2 A C13) 



(Cn + C12) C33 



2 (7^3 



(9) 
(10) 

Q 2 (11) 



7 (Cn - C 12 ) -2 /3 C_ 



(Cn — Ci 2 ) Cqg — 2 C 
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As (ei — e 2 ) and e§ display the same symmetry charac- 
teristic as Q, it is then natural to find that they are pro- 
portional to Q. We also obtain that e 4 = e$ = which is 
again compatible with a 4 — -> 2 symmetry change. Fur- 
thermore, using the minimization with respect to the or- 
der parameter Q, the corresponding expressions for the 
order parameter Q, the critical temperature T c and its 
pressure dependence can be written as 



(12) 
(13) 

(14) where 



Q(T,P) 
T c 

dP 



C a (T c + ^P - T) 



A 



C b 
a C a 



= 2 



C c 
a C a 



(15) 

(16) 
(17) 



C a 

c b 

C c 

A 



(Cn + C 2 ) C33 — 2 C 2 3 
7 2 (Cn - C12) + 2/3 09 C 66 - 2 7 Ci 6 ) 
S (Cn + C12 - 2 C13) + 2A (C 33 - C13) 
8SXC 13 - 2A 4 C 2 3 - 4A 2 C 33 + Cu(A 4 C 33 



26 2 ) + C 12 (A 4 C 33 - 25 2 ) 



(18) 
(19) 
(20) 
(21) 
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Finally, the elastic constants are calculated using 17 

8 2 F 



d 2 F fd 2 F\ 1 d 2 F 



de m de n dQde m \dQ 2 J de n dQ 
According to this model, the elastic tensor associated with the tetragonal 4 phase is 



(22) 



f Cxi 



A{T,P) 

r 4- P 

°12 + A(T,P) 
Cl3 




Pi 
A(T,P) 



Cl2 T A(T,P) 

r - ^ 

^ n A(T,P) 
Cl3 






\c 



16 



-Ci 



§2 
A(T,P) 



Cl3 








Cl3 








C33 











C44 











C44 












C16 - 
-C 



Pi 



A{T,P) 



16 1 yl(T,P) 





2 

/"( _ 7 
^66 A(T,P) 



(23) 



while the corresponding tensor for the monoclinic 2 phase is 



Cu 
C i2 4 

C13 - 



x O + A^c; 
2AQ^ Q 




7*Q + 
1"« 



C12 f 

C u 
C 13 + 

— C-IR 



X Q- 

2\QX Q 
W 







Cl3 " 

C 13 4 
C33 



Yq 
2\QX Q 

4,S 2 Q 2 



lX Q _ 
Yq 



Yq 




2~ t 5Q 
Yq 



C44 








+ 2 V Q 2 

CQ 








CQ 

C44 + 2t?Q 2 




C16 " 
— C16 



7^ Q+ \ 
, 7AfQ_ 



2-fSQ 
Yq 




c fifi -^ 



(24) 



with 



X Q+ 

Yq 



I3 + 2XQ 
/3-2AQ 

2A A Q 2 



J 2 (C 11 -C\ 2 )-2I3 2 C I 



66 



(C n 

r 



Ci2)C\ 



66 



2Ci6 



(25) 

(26) 
(27) 

(28) 



A rapid inspection of the elastic tensor, corresponding 
to the solution below T c (or above P c ), indicates that 
its form is consistent with the elastic tensor of the mon- 
oclinic class 2. This indicates that the coupling terms 
considered in the Gibbs energy includes all terms which 
are fully compatible with a 4 — » 2 symmetry change. 



ergy, the Landau energy, and the elastic energy as defined 
in Eq.0l- Eq.[S| and Eq.0 respectively. The main differ- 
ence being that the order parameter Q now corresponds 
to a strain combination defined as Q = a\(e.\ — e 2 )-\-Q,^e^. 
Thus, the coupling terms invariant with respect to the 
high temperature symmetry are 



C. Proper Ferroelastic Model 

In the case of a proper ferroelastic model, the deriva- 
tion is slightly different. We can still write the Gibbs en- 



FM) = SQ 2 e 3 + XQ 2 (ei + e 2 ) + ( e 4 e 5 Q + V [e 2 + e 2 ) Q 2 (29) 

Using the minimization of the Gibbs energy with respect to e», the solutions for e% + 62, &s, &i, &5, dT c /dP remain 
identical to those obtained for the pseudo-proper case while 

T =T -2C^ 6 + (Cn - Ci 2 )C 66 

c ° a(a 2 (C 11 -C 12 ) + 2a 1 (a 1 C 66 -2a 6 C 16 ) ' [ ' 
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Finally, as that the Gibbs energy is now only a function of strains, the elastic constants are directly obtained using 



Thus, in the case a 4 — > 2 proper ferroelastic transition, the elastic tensor for the paraelastic phase is given by 



/ Cn+alA(T,P) 
C±2-a\A{T,P) 


C 12 -ajA(T,P) 
Cu+a$A(T,P) 


C13 








C 16 + ai a 6 A(T,P) \ 


Cl 3 








-Cie - aia 6 A(T, P) 


Cl3 


Cl3 


C33 




















C44 




















C44 





V Ci 6 + a ia6 A(T, P) 


-Ci 6 - aia 6 A(T, P) 











C 66 + alA(T,P) / 



(32) 



where A(T, P) = a(T-T c - %±P). Thus, the temper- 
ature and pressure dependence of the elastic constants 
should be linear. This differs from the pseudo-proper 
model where the variation in the elastic constants is ex- 
pected to be inversely proportional to A(T,P). Nat- 
urally, depending on the strength of the coupling pa- 
rameters this dependence could be quasi-linear. Thus, 
a clearly distinction between these two models can only 
be obtained through a detailed comparison with experi- 
mental results. 
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